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We study the interplay of superfluidity and glassy ordering of hard core bosons with random, frus- 
trating interactions. This is motivated by bosonic systems such as amorphous supersolid, disordered 
superconductors with preformed pairs and helium in porous media. We analyze the fully connected 
mean field version of this problem, which exhibits three low temperature phases, separated by two 
continuous phase transitions: an insulating, glassy phase with an amorphous frozen density pattern, 
a non-glassy superfluid phase and an intermediate phase, in which both types of order coexist. We 
elucidate the nature of the phase transitions, highlighting in particular the role of glassy correlations 
across the superfluid-insulator transition. The latter suppress superfluidity down to T = 0, due to 
the depletion of the low energy density of states, unlike in the standard BCS scenario. Further, we 
investigate the properties of the coexistence (superglass) phase. We flnd anticorrelations between 
the local order parameters and a non-monotonous superfluid order parameter as a function of T. 
The latter arises due to the weakening of the glassy correlation gap with increasing temperature. 
Implications of the mean fleld phenomenology for finite dimensional bosonic glasses with frustrating 
Coulomb interactions are discussed. 



PACS numbers: 

I. INTRODUCTION 

Bosons in random environments occur in a variety 
of experimentally relevant systems, ranging from cold 
atomic gases, superconductors and quantum liquids. The 
superfluidity of Helium-4 ('*He) in porous media was one 
of t he f irst phenomena observed in this type of sys- 
tems,'il2l featuring an interesting competition between 
Bose-Einstein condensation and localization by random 
potentials.^ 

In more recent years, supersolidity in crystalline '^H^^ 
has been reported. It soon became clear that defects of 
the crystalline order and amorphous solids sustain more 
robust supersolidity, spurring the idea that disorder, or 
even glassy order, may be a crucial element in under- 
standing the superfluid part of those systems.'^Hli' 

A recent experimenlP^ reported indeed that the super- 
solidity in *He is accompanied by the onset of very slow 
glassy relaxation. This suggested that an amorphous 
glass with a superfluid component is forming, a state of 
matter which was dubbed a "superglass". These exper- 
imental results have motivated several theoretical inves- 
tigations into the possibility and nature of such amor- 
hously ordered and yet supersolid systems.^i^lti^ 

Similar questions as to the coexistence and interplay of 
glassy density ordering and superfluidity arises in disor- 
dered superconducting films, which feature disorder- or 
field-driven superconductor-to-insulator quantum phase 
transitions.^^ In several experimental materials, this 
transition appears to be driven by phase fiuctuations 
of the order parameter rather than by the depairing 
of electrons, suggesting that the transition can be de- 
scribed in terms of bosonic degrees of freedom onlyf^ 
This had lead to the dirty boson modeP^^ and the no- 



tion of the Bose glass,!^ in which disorder and interac- 
tions lead to the localization of t he bosons, while the 
system remains compressible i^ this context, the 
term "glass" refers mostly to the amorphous nature of 
the state rather than to the presence of slow relaxation 
and out-of equilibrium phenomena due to frustrated in- 
teractions. However, if superconductivity develops in a 
highly disordered environment, such frustration may add 
in the form of Coulomb interactions between the charged 
carriers, which may become important, as screening is 
not very effective. It is well known that in more in- 
sulating regimes, strong disorder and Coulomb interac- 
tions ma y in duce a glassy state of electrons (the Coulomb 
glass). l^SHSZl It is therefore an interesting question whether 
such glassy effects can persist within the superconduct- 
ing state of disordered films. A memory dip in resistance 
versus gate voltage (similar to the conductance dip in in- 
sulators, which is considered a smoking gun of electron 
glassiness^^) was reported as a possible indication of such 
a super-glassy state in Ref. UHl even though doubts about 
its intrinsic nature were raised later on.l^Sl 

The recent developments in ultra cold atom^^D open 
new ways to studying bosonic atoms in the presence of 
both interactions and disorder. Those can exhibit su- 
perfluid or localized, and potentially also glassy phases, 
especially if the interactions are sufficiently long ranged 
and frustrated, as is possibly the case for dipolar inter- 
actions. 

Motivated by these experimental systems, we study a 
solvable model of bosons subject to disorder and frus- 
trating interactions, as proposed previously in Ref. 1151 
This solvable case provides insight into the possibility of 
coexistence of superfiuidity and glassy density order, as 
well as into the nature of the coexistence phase (the su- 
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perglass). In particular, for the considered mean field 
model we prove the existence of a superglass phase. This 
complements the numerical evidence for such phases pro- 
vided by quantum Monte Carlo investigations in finite 
dimensions'^ and on random graphs.'^ Those were, how- 
ever, limited to finite temperature, and could thus not 
elucidate the structure of the phases at T = 0. In con- 
trast, the analytical approach allows us to understand 
the quantum phase transition between glassy superfluid 
and insulator, and the non-trivial role played by glassy 
correlations. Indeed it is an interesting feature of the 
considered model that despite of its infinite connectiv- 
ity, it features an insulating, Anderson localized phase 
at sufficiently small but finite hopping strength. This is 
in contrast to non-glassy models where the localization 
and insulating behavior is lost in the large connectivity 
limit, unless the hopping is down-scaled logarithmically 
with the connectivity.'SS 

The model studied here is an interesting type of a 
quantum glass.l^ Like other canonical quantum glass 
models of mean field type, such as the random exchange 
Heis enber g model, the Ising spin glass in transverse 
ggy|35]36l g^jjj frustrated rotor models,'^^!^! the glass 
phase of the superglass model breaks ergodicity, and will 
exhibit a large number of metastable states and associ- 
ated slow relaxation and out-of-equilibrium dynamics. It 
has recently been pointed out that such mean field glasses 
can faithfully be realized by atoms in laser cavities, which 
are coupled at long distances through interactions via dis- 
crete photon modes. "^^ *' 

In the above mentioned canonical models the non- 
glassy phase is usually quantum disordered, having no 
broken symmetry. However, a new situation arises in 
the model studied here: The glassy ordering competes 
with a different order: the superfluidity, or the transverse 
(XY) ferromagnetism in a magnetic analogon. The mi- 
croscopic coexistence of these two types of order, which 
we demonstrate to exist in our model, is rather non- 
trivial, as one might instead expect a first order transi- 
tion between phases with either of the two orders. At 
zero temperature, the quantum glass transition is in- 
duced by frustrated interactions, which win over weak 
quantum fluctuations. While the quantum phase transi- 
tion of such models is understood relatively well,^^ the 
deep glass phase wit h its c ollective gapless modes has not 
been much explored.'^^'^ll 

Let us finally remark that the glassy, amorphous super- 
solid, which the "superglass" phase constitutes, is quite 
different from the type of supersolid proposed theoreti- 
cally in the early seventies.!^ In those scenario the bosons 
organize spontaneously on a lattice, which breaks trans- 
lational symmetry, but is incommensurate with the bo- 
son density, allowing for vacancies to move through the 
solid. Our model considers instead bosons on a prede- 
fined lattice, on which an inhomogeneous density pattern 
establishes in the glass phase. 

The remainder of the paper is organized as follows: In 
Sec. II we introduce the mean field model of a superglass. 



We obtain the effective action for a single spin with both 
the replica method and the cavity approach in Sec. Ill, 
and introduce the key concept of the distribution of local 
fields. The self-consistency equations of the mean field 
theory are solved under static approximation, which is 
argued to be exact in parts of the phase diagram. Sec. IV 
determines the instabilities towards forming glassy and 
superfluid order, and establishes the phase diagram fea- 
turing two quantum phase transitions at T = 0. We also 
discuss the robustness of the phase diagram to random 
potential disorder. In Sec.V, we study the bulk of the 
"superglass" phase in more detail. We show that the 
glass and superfluid order parameters are locally anti- 
correlated. Moreover, we flnd the superfluid order pa- 
rameter to have an interesting non-monotonous behavior 
as a function of T. The implications of this mean field 
analysis for realistic finite dimensional models, e.g. with 
frustrating long range Coulomb interactions, will be dis- 
cussed in Sec. VI. Some detailed derivations are collected 
in the appendix. 

II. MODEL 

We consider the fully connected model of hard-core 
bosons with random pairwise interactions between all 
bosons, 

= - E + E ^^"^ ~ ^ E (^l^j + ^h) •(!) 

i<j i i<j 

Here, ni = b\bi is the number operator on site i, and the 
hard core constraint limits Ui to assume values or 1. 
bi{b\) denote the annihilation (creation) operators for a 
hard-core boson at site i. Vij is a quenched disorder with 
Gaussian distribution of zero mean and variance V"^ /N , 
ti describes a quenched disorder potential for the bosons, 
and th/N is the unfrustrated hopping strength between 
any pair of sites. The scaling of the couplings with is 
chosen so as to yield a non-trivial thermodynamic limit 
for A^ — > oo. 

In the absence of hopping, the model becomes classi- 
cal and is equivalent to the Shcrrington-Kirkpatrick spin 
glass (SK) model^® in a random field. It is well-known 
that except for lowering the transition temperature, the 
random fields do not alter the low temperature proper- 
ties of the glass phase. We thus restrict our attention 
mostly to a slightly simpler model proposed in Ref. I15[ 
which corresponds to a special choice of the e^: 

H = -Y.V.,,{n,-l/2){n,-ll2) 
i<j 

-^E(^I^.+^.S- (2) 

i<j 

Similar fermionic mean fleld models have been studied 
in Refs. I49I50I The identiflcation 2n, - 1 = sf e {±1}, 
bl = sf ,bi = s~ , allows us to map this model into a fully 
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connected spin glass model with quantum fluctuations 
arising from non-random spin flip terms, 



with 



H 



1<J 



N ^ 



(3) 



t = 



(4) 



with the simple dictionary 

For i = 0, this Hamiltonian reduces to the SK model, 
which possesses a spin glass phase at low temperature, 
T < Tg = J. Without the Ising interactions, Jij — 0, 
the Hamiltonian turns into the mean field XY model, 
which has a superfluid (or XY ferromagnetic) phase at 
low temperatures [T < Ts = t). In this paper we estab- 
lish the phase diagram and study the properties of the 
bulk phases resulting from the competition of random 
density-density interactions and boson hopping (bosonic 
language) or equivalently random Ising interactions and 
ferromagnetic transverse coupling (spin language). 



III. FREE ENERGY AND SELF-CONSISTENT 
EQUATIONS 

A. General formalism 



The disorder average of the free energy of the model 
(Isl) can be obtained using the replica methocP^ 



(logZ)^ 



lim 



(5) 

indicates an 



where Z is the partition function and (. 
average over the couplings Jij. 

Following a method introduced by Bray and Moor^^ll 
it is useful to represent the partition function as an imag- 
inary time path integral: 



Z" Trrexp<^/3 




'0 a=l i<j 



(6) 



where T orders the operators in decreasing order of their 
argument r € [0, 1]. This "time" argument of s(t) merely 
serves us to define the time-ordering, while s(t) denotes 
always the same Pauli matrix, independently of time. 

Averaging over disorder and decoupling the spins on 
different sites using a Hubbard-Stratonovich transforma- 
tion with the order parameter fields Qab, , M^, we 
obtain: 



l[dQaa{T,T')dM:{T)dMy{r) 



xl[dQab{T,T')exp{-Nr) (7) 

a<b 



2o2 rl rl 



"'0 



drdr' 



2 



logZ, 



M^irr + Myirr 



Z = TrTcxp (— S'cft) ■ 



ScS — — 



^0 



drdr' 



(8) 



(9) 



(10) 



-tPY.[ dr (^M:(r).^(r) + Afi'(r).^(r)^ . 

In the limit iV — > oo, the functional integral ([T]) is dom- 
inated by the saddle point of the replicated free energy 
J^, which satisfies 







5T 



SQabir, t') 

{sl4)cS = Qab 
5Qaa{T,T') 

R{t,t'), 



Qab{T,T')^{Tsl{T)sl{r')}^^ 



(11) 



,(r,r') = (rs:(r)s:(r')),ff 



(12) 



(r) = Kir))^^ ^ M\ (13) 



^ My{T) = (slir))^^ ^ My, (14) 



SJ' 



where denotes the average with respect to the ef- 

fective action ScS of a single site. We have used that, 
as usual, the saddle point values of Qa^b and Mf'*' 
are independent of imaginary time, while QaaiT, r') de- 
pends only on the imaginary time difference.SZl Further- 
more, Qaa and Mf^'y do not depend on the replica index 
a. For Qab we make the standard ultrametric ansatz, 
parametrized by a monotonous function q{x) on the in- 
terval X & [0, l]j which is well-known to describe success- 
fully the SK model and other mean field glasses.HSWe are 
free to choose coordinates in the x, y plane such that the 
spontaneous magnetization M points in the x-direction, 
and thus we set = 0. 

Note that ^ signals the presence of transverse 
(XY) order of the spins, that is, superfluidity of the hard 
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core bosons, which breaks the U{1) symmetry sponta- 
neously. On the other hand, a non-constant value of Qa^b 
implies the spontaneous breaking of the replica symme- 
try, and thus the presence of a glass phase with many 
metastable states and non-trivially broken ergodicity. As 
long as we do not consider random field disorder, the 
breaking of replica symmetry coincides with the breaking 
of the Ising symmetry and is signalled by a nonzero value 
of Qa^b- We will see below that the t/(l) and the replica 
symmetries can be broken simultaneously in a what has 
been called a "superglass phase" in Refs. 1151151 

To find the location of a (continuous) glass transition, 
we expand the free energy to second order in Qat- We 
find an instability towards replica symmetry breaking, 
and thus the emergence of a glassy density ordering of 
bosons, when 

13 J C t drdr' {Tsl{T)sl{r'))^^ = /3J C drR{T) 



^0 







1, 



or 



(15) 



(16) 



where x = Xzz{^ = 0) is the zero- frequency limit of the 
longitudinal susceptibility. This condition is of course to 
be evaluated at Qa^b — 0. 

On the other hand, a second order phase transition 
from the high temperature phase towards a superfluid 
state is indicated by the instability condition, which fol- 
lows from d'^F/dAP = 0: 



^0 



drdr' {Tsl{r)sl{T')) 



off 



(17) 



B. Solution of the saddle point equations 

A full solution of the saddle point equations involves 
the solution of the problem of interacting replica as well 
as the evaluation of dynamical correlation functions with 
the effective action S'efi, if M ^ and the replica sym- 
metry is broken as well. 

Here we describe what steps an exact solution involves, 
and then discuss the approximations we will use to study 
parts of the phase diagram, especially the bulk of the 
superglass phase. 

To describe a non-glassy superfluid phase, the replica 
structure is trivial, and one needs to solve the self- 
consistency equations 



M 



/cff ' 



R{T)^{Ts^{T)s^m.^: 



(20) 



with effective action 



B'^ P 

SoS = / / dTdT's'{T)R{T~T')s'{T') 



'0 JQ 

-fitM [ dTs'^iT). 
Jo 



(21) 



These can be solved using techniques as used in dynam- 
ical mean field theory.'^^ 

In a glassy phase the replica structure has to be taken 
into account. Assuming the standard ultrametric struc- 
ture of the saddle-point matrix Qab, the above single- 
replica scheme has to be generalized to include a self- 
consistent distribution of frozen longitudinal fields P{y) 
acting on a given replica. This captures the distribution 
of random frozen fields yi created by the exchange of 
sites i with the frozen magnetization pattern with a spin 
glass state.!^ In practice this requires the simultaneous 
solution of 



or 



tX^ - 1, 



(18) 



where x'^ = Xxx{^ = 0) is the static transverse suscepti- 
bility. These expressions must be calculated in the non- 
superfluid phase where M = 0. In this regime the effec- 
tive action Scs is classical, which entails the further sim- 
plification R{t) = 1. This feature is due to the suppres- 
sion of quantum fluctuations in the non-superfluid phase 
by factors of 1/A'', due to the scaling of the transverse 
coupling. It allows us to find the superfluid-insulator 
transition analytically, even at zero temperature, with- 
out solving a full quantum impurity problem. In partic- 
ular, we immediately find that the transition from the 
disordered high temperature phase to a glassy phase is 
given by 



(19) 



exactly as in the classical SK model. However, the glass 
transition line will be modified if it is preceded by a su- 
perfluid transition at higher temperature. 



■m{y) = 




■nixiy) = 




M = 


j dyP{y)mx{y), 


R{t) = 


J dyP{y) {Ts^T)s^{Q))s^^^^y^ 



where the effective single replica action in a frozen field 
y reads 



Scs{y) = -^^-^ [ [ drdr' 
^ Jo Jo 



XS-(t) [R{T-T')-qEA]s^T') 

-PtM ( dTs'^ir)- (5y [ dr s'- (t) . {23) 
Jo Jo 

The Edwards- Anderson order parameter ^ea character- 
izes the glassy freezing in a pure state of the glass and is 
given by 



gEA = g(a; = 1) = / dyP{y)m'{y). (24) 
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As first derived by Somniers and Dupontf^ the frozen 
field distribution P{y) = P{y,x = 1) is obtained from 
a self-consistent solution of the differential equations on 
the interval x E [0,1] 



'm{y,x) = - 
P{y,x 



q(x) 
2 

q{x) 



m"{y, x) + 2xm(y, x)m'{y, x) 
P"{y, x) - 2x{m{y, x)P{y, x))' 



, (25) 
,(26) 



with 



i{y, x^l) = m{y), P{y, x = 0) = S{y), (27) 



where dots and primes denote derivatives with respect 
to X and y, respectively. The solutions of these differ- 
ential equations solve the saddle point equations of the 
replica free energy ([s]).'^ The overlap function q{x), which 
parametrizes the ultramctric matrix Qab by the distance 
X between replica, must obey the self-consistency relation 



q{x) = / P{y,x)m{y,xf 



(28) 



Notice that Eqs. ( [25| ) and (26 1 are the same as in a clas- 
sical spin glass. The influence of quantum fluctuations 
enters through the boundary condition m{y,x = 1) = 
'm{y), where m{y) was defined in Eq. (22) These differ- 



ential equations provide an elegant way of integrating out 
all spins except for one.*^ 

Once the above scheme has been solved self- 
consistently, site-averaged observables such as the lon- 
gitudinal magnetization are given by 



dyP{y) i^lsMv) 

dyP{y)m{y). (29) 



The properties of the solution of these differential equa- 
tions are well understood in several classical models ex- 
hibiting full replica symmetry breaking with continuous 
functions q(x)¥^^^Th.e full solution of mean field quan- 
tum glasses in the ergodicity broken has not been ana- 
lyzed in the literature so far. However, an analysis of 
the transverse field SK model shows that most features 
of the low temperature solution of q{x) carry over rather 
naturally to the quantum case.l^ A salient new feature 
in the quantum case is the fact that full replica sym- 
metry breaking implies marginal stability of the whole 
glass phase, which in turn ensures the presence of gap- 
less collective excitations. The latter is very similar to 
what was found, e.g., in the threshold states of quantum 
p-spin models,'^ or in the quantum dynamics of elastic 
manifolds, approximated with a replica symmetry break- 
ing variational approach.®*^ 



C. Alternative derivation by a cavity approach 

The replica-diagonal part of the above scheme will be- 
come easier to understand, if we derive it in a cavity 
frameworlPSl similarly to the derivation of the quantum 
analog of Thouless- Anderson-Palmer equations by Biroli 
and Cugliandolo.^^ From a cumulant expansion in the 
couplings involving site o it is easy to obtain the follow- 
ing effective action for the site o: 



QOff _ 



(3^ 



2 



JQ 



drdr' 



sl{r') 



-/3 dT\hl{T)sl{T) + h^{r)sl{T)\. (30) 



Here 



'Xr)r = Y,jo^ {si 



(31) 



is the site-dependent longitudinal field, which does not 
depend on time, however. The index o denotes a "cavity 
average", i.e. an average over the action of the system, 
in which the site o has been removed. The subscript c 
indicates a connected correlator. The effective transverse 
field, 



N 
t 

N 



E 



sf ) = tM, 



(32) 



does not fiuctuate from site to site, and is independent of 
T if we neglect subleading terms, which scale as inverse 
powers of N . Note that for large N 



<1EA\ 



(33) 



independently of the site o. The distribution of hf over 
the sites i is the frozen field distribution. 



N 



Piy)^N-'J2^iy 



(34) 



i=l 



computed in the replica formalism. Thus we precisely 
recover the self-consistency problem for the replica diag- 
onal, while the solution of the replica off-diagonal part 
furnishes the distribution P{y). 

For the study of the phase transition from the insulat- 
ing glass phase into the superfluid, it will prove crucial 
to use the full low temperature solution of the SK model. 
However, in order to analyze properties of the mixed "su- 
perglass" phase we will restrict ourselves to a one-step 
approximation, which we discuss in the next section. 
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D. Static and one-step approximation 

In order to avoid solving numerically a full self- 
consistent quantum problem as outlined in Eqs. ( 22 ) 



above, we will resort to the widely used static approxima- 
tion. The latter consists in seeking a minimum of the free 
energy not with respect to the full function space R{t) 
but, instead with respect to a constant value R{t) — > R. 

A further approximation which we will use in the study 
of the quantum glassy phases is the one-step approxima- 
tion for the structure of replica symmetry breaking. It is 
equivalent to assuming a step form of q{x) 

q{x) ^ e{x - xi)Qi + Qixi ~ x)Qo (35) 

and optimizing the free energy over xi,Qi,Qo. This is 
expected to give qualitatively good results, especially at 
intermediate temperatures and close to the glass tran- 
sition. Combined with the static approximation for the 
replica diagonal, short time part one obtains the free en- 
ergy functional per spin: 

4 

^ - Dyo 



[ix^^l)Ql~xiQl+R^] 



(36) 



2 



Xl 



X log / Dyi 



DyR2cosh{/3^hl+t^M^) 



where h,,, 



yo + yi + VR- Dyo, Dyi and Dya are 
f yl 

exp ( - ^^-^ 

Gaussian measures: ZJj/o = 



iQa-f 



^dyo, Dyi 



dyi and Dya 



dyn- 



Note that Qi — (?ea is the Edwards Anderson order 
parameter in the one step approximation, while Qq is 
the overlap between different spin glass states. We point 
out that the above free energy differs from the expression 
given in Ref. 1151 where the static approximation was not 
carried out correctly. This error was at the origin of 
several strange features of the phase diagram reported 
there, such as a T-independent transition between su- 
perfluid and superglass and a J-independent superfluid 
transition. 



In the case of one-step replica symmetry breaking, the 
frozen field distribution within one pure state can be 
obtained by stepwise integration of the flow equations 
( p5p6| , yielding (cf. 



Piy) - I i^.o ^^^^^^^,:LV^f^^V^^i^°^^^\ (38) 

J j DyiZ^^iyo + yi) 

where Djji is a Gaussian measure like Dyi with variance 

(Qi-Qo)J'. 



Under the static approximation, Eq. (37) becomes 



where 



ZM-jDynZ.Uy + yn), 



,{y) ^ 2 cosHfi ^ y^ + MH^). 



(39) 



(40) 



One can interpret yn as a random field, which is gen- 
erated by the thermal fluctuations of the non-frozen part 
of the magnetization. 

The longitudinal and transverse magnetizations of a 
spin in a frozen field y introduced in Eqs. (22) are easily 
seen to be given by 



m{y) 



-J-\og{Z,^{y)), (41) 
/? dy 



The saddle point equations for the Edwards- Anderson 
parameter Qi and the superfluid order parameter M can 
now be expressed as: 



dyP{y)vr?{y), 



(43) 



= ^ E = / dyP{y)m.{y). (44) 

i 

The saddle point equation for the parameter R reads 



E. IRSB free energy and self consistent equations 

Here we rewrite the one-step self-consistency equations 
with the help of the local field distribution. 

The effective partition function of a single spin is 

Z,s{.y) = Trrexp(-5eff(y)) (37) 
TrTexp ( I I drdr' 



'Q JQ 

s^r) [R{T-T')-qEA]s^T') 
+l3tM [ dTs''{T) + l3y [ dTs'ir) 



m-Qi)= / P{yWL{y)dy^ 



(45) 



which relates the static approximation of the connected 
s^-correlator, i? — Qi, to the average local susceptibility 



(I , . 5m(y) 



(46) 



The saddle point equation for the Qq can be written 
in a similar way: 



Qo^ dyoP{yo;xi)m'^{yo;xi), 



(47) 
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where 



m{yo;xi) 



exp( 



yi)m{yQ 



yi) 



J DyiZ^^{yo + yi) 



(48) 



are discrete versions of the continuous functions 
P{y,x),m{y,x) introduced above. 

Optimizing the one-step free energy with respect to 
Qi, M, R and Qq yields the saddle point equations 
Eqs. ( 43]|47 ). To capture equilibrium states, we should 
further extremize with respect to the Parisi parameter 



xi, I.e. 



df. ^ 
dxi 



0, which yields the further condition 



xi / Dyo 



(Ql ~ QlW = J Dyolog J Dy,Z^^{yo + yi) 
J DyiZ^^{yo + yi)\og Zcffiyo + yi] 



jDy.Z^^iyo + yi) 



(49) 



It is a useful check that upon imposing Qi = Qo, the 
saddle point equations for Qo and Qi reduce to the same 
replica symmetric constraint. When M — 0, the local 
field distribution, the free energy and the saddle point 
equations reduce to those of the classical SK model, as it 
should be. 



IV. PHASE DIAGRAM 

Let us now study the phase diagram of our model 
The gross features of the phase diagram we find are sim- 
ilar to the ones found in Refs. Il5fl61 The low tempera- 



ture phase exhibits three phases: a non-glassy supcrfluid 
at small J/t, an insulating (non-super-fiuid) glass phase 
at large J/t, and most interestingly, a phase in between 
with both glassy order and superfiuidity. However, as 
mentioned before, we find a distinctly different behavior 
of the phase boundaries than Ref. [151 

Moreover, we are able to analyze the limit T — » 0, 
whose properties were inaccessible in previous works-^^^E^ 
The latter is of particular interest in the context of the 
superfluid-insulator transition. 

The findings of the mean field analysis are in quali- 
tative agreement with Monte Carlo studies in finite di- 
mensions at low but finite temperatures. The analytical 
approach allows for a detailed analysis of the properties 
of the mixed phase, and of the glass-to-superglass tran- 
sition. 



A. High temperature phase 

The high temperature phase is simple to describe. 
Since M — 0, the system behaves identically to the para- 
magnetic phase of the classical SK model, and R = I 
holds exactly. In this regime the static approximation is 
of course exact. 




Fig. 1: Phase diagram of glassy hard core bosons. At high 
temperature, the straight blue line T — J indicates the clas- 
sical SK glass transition line. The red solid line shows the 
superfluid phase boundary which is given by the instability 
condition (51 1. The two lines cross at the tricritical point 
{T/t)T = {J7t)T = 0.7248. At low temperature, the blue hue 
shows the phase boundary of the glass within the superfluid 
phase, as evaluated within the static approximation, cf. ( |16[ l . 
The glass transition at T = occurs at (J/t)g, stat = 1/2. 
The red solid line indicates the location of the onset of su- 
perfluidity within the glass phase, as evaluated within the 
full breaking of the replica symmetry to the instability con- 
dition ( |63| |. The superfluid transition at T = takes place at 
{J/t)s = 1.00. 



At large enough J/i, the leading instability upon low- 
ering the temperature is the classical glass transition at 
Tg — J, as mentioned earlier. However, at small val- 
ues of J/t, the tendency to form a superfluid wins. The 
instability condition towards XY symmetry breaking. 



dm^{y = 0) 



dhx 



1 



(50) 



can be evaluated exactly. In this expression hx is a 
uniform transverse field. The transverse susceptibility 
is easily calculated for the replicated Hamiltonian with 
a Hubbard-Stratonovich transformation of the quadratic 
term R J drdr' s'^{t)s^{t'). This results in the instability 
criterion 



1_ ^ _ Jdhe-^ /^-^ sinh(/?fe)/^/i 
~^ /d/ie-^V2J^cosh/3/i ■ 



(51) 



The glass transition and the superfluid transition line 
cross at the tricritical point 



Tt — Jt — tx 



f dze ^ sinh(z)/z 
/ dze'~^'^ /"^ cosh z 



0.7248 iT- (52) 



The result does not have the familiar looking form 
of an average local transverse susceptibility. However, it 
can indeed be recast in such a way. This furnishes us 
a better understanding of the interaction effects in the 
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high temperature phase, and at the same time illuminates 
the nature of the static approximation in the superfluid 
phases. 

Let us rederive the above result directly from the non- 
replicated Hamiltonian: 



(53) 



where is an infinitesimal field. In a given classi- 
cal Ising configuration, the spin i sees an "instanta- 
neous" local field hi — J2j^iJij^j' while the trans- 
verse coupling is negligible in the paramagnetic phase 
where N-'^J2j{sj) = M = = ^"^Ej(si)- Thus 
the transverse susceptibility can be calculated as a site 
and configuration average of the susceptibility of a sin- 
gle spin sitting in an instantaneous field h, X^(/i) — 
/p^dr(s^(T)s=^(0)) ^t&nh{l3h)/h. 

The thermal distribution of instantaneous local fields 
of the SK model has been well studied,'^ and takes the 
rather simple form 



Pinst(/l) =C0sh(/3/l) 



exp(- 



'2.P 



2 



) 



(54) 



in the paramagnetic phase. Note that the instantaneous 
field distribution is not a simple Gaussian, but small 
fields are under-represented. This phenomenon is closely 
related to the suppression of small fields encountered in 
the cavity approach to Ising systems, and is a precursor 
effect of the openin g of the pseudogap in glassy phases 
at low temperaturesPSMMl 

The total transverse susceptibility is obtained as an 
average of the local susceptibihty x(^) over Pinst{h): 



dh 



s'mh l3h 



\/27rJ2 



(55) 



which indeed coincides with the replica result (51 ) 



The static approximation for superfluid phases has a 
completely analogous effect. The approximation replaces 
the dynamically fluctuating exchange fields on the vari- 
ous sites by a random distribution of quasi static fields. 
The latter differs from the distribution of frozen fields 
(which is Gaussian at high T) by a random Gaussian 
smearing with variance J2(i? — qea), and a reweighing 
factor proportional to cosh(/3/i) which accounts for the 
fact that a small instantaneous fields is less likely to ob- 
serve on a given site, as it implies a positive free energy 
fluctuation in the environment. 



B. Onset of glassy order within the superfluid 

The instability towards forming a glass occurs when 
Jx^l = I3JJ dTdr'RiT - r') = 1. Within the super- 
fluid phase it is difficult to calculate this susceptibility 
exactly, and we thus first resort to the static approxima- 
tion, R{t — t') — > R. The instability of the statically ap- 
proximated free energy occurs when J PR = 1, PR being 
the static approximation for the longitudinal susceptibil- 
ity x" . Within the non-glassy superfiuid phase there are 
no frozen fields, P{y) = 5{y). Thus, from Eqs. ( [44p5| , 
the two relevant saddle point equations read 



M = m^(y = 0), 



(56) 



(57) 



where m^ly = 0) and xJLc(j' ~ ^) ^""^ evaluated 
from Eqs. ( [39p2l ) and Q. 

They have a relatively simple low temperature limit. 
One verifies that it is self-consistent to assume that 



PR = x" 



t 



M 



1 — ra- 



T 
i' 



(58) 



with finite numbers r, to, as T — > 0. 

Injecting this into the above self-consistency equations, 
and evaluating the Gaussian integral over yn in Eq. ( 39 ) 
around the stationary point, the equations simplify to: 



= 1 



t2 - J2r 
J2r 



2(i2 - r J2) 



f o(rA), 

0{T/t). 



(59) 
(60) 



This yields the solution for the susceptibility Jx" 



l-4(f) 



The static approximation predicts the quantum glass 
transition to occur at the critical point 



(;,stat 



1 

2' 



(T = 0) 



(61) 



where Jx" = 1. 

It is difficult to predict whether we over- or underesti- 
mate the phase boundary with the static approximation 
in the superfiuid phase. This is because the approxima- 
tion has two competing effects with respect to the onset of 
glassy order. On one hand, we approximate the dynamic 
longitudinal susceptibility by the static one. Since the 
latter is bigger, we tend to overestimate the stability of 
the glassy ordering of . This effect is well-known from 
the SK model in a (constant) transverse field r.EzMmiOn 
the other hand, the static approximation underestimates 
quantum fluctuations of Sx, at least at low T . Indeed we 
see above that at T = 0, the static approximation pre- 
dicts maximal transverse order, M = 1, independently of 
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the vahie of J/t, while it is easy to show that quantum 
fluctuations around the transverse ferromagnetic state 
decrease the magnetization as A/ = 1 — 0{{J/t)'^). The 
overestimate of M leads to an underestimate of the longi- 
tudinal susceptibility, and thus of the tendency to glassy 
order. In view of these competing tendencies, it is hard to 
predict on which side with respect to Eq. (61 1 the exact 
glass instability will be located. 

However, there is a simple way to obtain an upper 
bound for the quantum critical point. In the superfluid 
phase our model ([s]) is very similar to the SK model in 
a constant transverse field with the difference that 
the effective transverse field Mt is self-generated and has 
to be determined self-consistently. However, it is clear 
that the effective transverse field is always smaller than 
t. From quantum Monte Carlo results for the transverse 
field SK model, one knows that a quantum glass phase 
obtains for J/T > 0.76.1^ This implies that the model 
studied in the present work must certainly be in a glassy 
phase if J/t > 0.76. The latter value is thus an upper 
bound for {J/t)g. Approaching from large values of J/T 
we will find below in Eq. ( 66 ) that the non-superfluid 



field distribution as 



glass phase becomes unstable towards superfluidity al- 
ready at {J/t)s = 1-00. Hence, we conclude that a phase 
with both superfluid and glassy order parameters exists 
for a substantial range of parameters covering at least 
the interval 0.76 < J/t < 1.00. 



C. Superfluid instability within the insulating glass 

phase 

1. Instability criterion 

Our discussion of the phase boundaries will be com- 
plete, once we have addressed the superfluid instability 
with in the glass phase at large J/t. The instability con- 
dition reads 



t / dyP{y) 



1, 



(62) 



where P{y) is the non-trivial distribution of frozen lo- 
cal fields in the classical glass phase of the SK model. 
The properties of P{y) are well studied, and turn out 
to be crucial to understand the low temperature behav- 
ior of the phase boundary and the physics of the glassy 
superfluid-to-insulator quantum phase transition. 

We recall that in the non-superfluid glass phase the 
static approximation is exact with R ~ 1, so that the 
instability criterion can be expressed in the form. 



^ J dyP{y) 



= 1, 



(63) 



DyRsm\v{P{y + yR))^ 
jDyn cosh(/3(y 2/ii)) 

exp(- 27^3^1^). This con 



where Dyo = , i^^p, — 757^ ^ 

^" ^27r(l-gEA)J2 2(l-qEA).. 

dition can be expressed in terms of the instantaneous 



t J dhPinst{h) 



tanh(^/i) 
h 



where 



Pi„st(/i) = / Piy)dy 



I 



cosh /3 ft, 



exp I 



2ho 



I3hc 
1 



cosh I3y ^2TTho/l3 



(64) 



(65) 



is the instantaneous field distribution, which was first de- 
rived in Ref. 61. The term ho — /3J^(1 — (7ea) is known as 



Onsager's back reaction. Eq. (64 1 can be recognized as a 



BCS-equation, where the instantaneous field distribution 
Pii\st{h) takes the role of the density of states. 

The temperature dependent local field distribution 
can be obtained from a numerical solution of the self- 
consistent set of full RSB equations ( 25][28 1, from which 
the phase boundary of the insulator-to-superfluid transi- 
tion is deduced . This yields the solid [red] line in Fig. [I]). 
For comparison we also evaluate the phase boundary 
within a one-step approximation, which works well at 
moderate temperatures. However, it fails badly at low T 
where a non-physical reentrance of the superfluid insta- 
bility would be predicted, and the quantum phase tran- 
sition at T — is completely missed. 

We note in passing that the thermodynamics of the in- 
sulating phase is essentially classical because of the scal- 
ing of the transverse coupling as t/N . If instead t were 
random and scaled as 1/ VN , the glass phase would also 
exhibit quantum fluctuations and would not reduce to 
the purely classical SK model. In that case, the analysis 
of the transition would become much more complicated. 
However, even though the thermodynamics can be ob- 
tained by a purely classical saddle point computation, 
one should not conclude that excitations do not have any 
quantum dynamics. 



2. Low temperatures and quantum phase transition 

At low temperatures, the most prominent feature of 
the local field distribution P{y) is a linear pseudogap 
which opens at small fields. The latte r is r equired to 

in a very 



I64I65I 



assure the stability of the glass phase 
similar manner as the Efros-Shklovskii Coulomb gap 
arises in electro n gla sses with unscreened, long range 
1/r interactions.'^^'^ More precisely, it is known that 
P{y) ^ a\y\ + 0{T) with a = 0.301 for fields in the 
range T <C |y| <C J , while the distribution decays like a 
Gaussian for \y\ ^ J. At zero temperature the pseu- 
dogap extends down to y = (i.e., the chemical po- 
tential in the terminology of hardcore bosons), while at 
finite but low temperatures T <^ J, P{y) assumes a 
scaling form P{y) = Tp (y/T ) with P(0) — const, and 
p{x > 1) = a\x\ + const. .E^EZl 

This scaling form asserts that only a fraction of {T/J)"^ 
is thermally active. Therefore the Edwards- Anderson pa- 
rameter tends to 1 as 1 — q-^A ~ [T / J)"^ ■ Accordingly, 
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as r — there is no difference between the distribu- 
tion of frozen and instantaneous fields, Pinst(^), since no 
thermal fluctuations are left. In this limit the instability 



condition (63) for onset of superfluidity then takes the 
form: 



dy 



\y\ 



1. 



(66) 



Using the above mentioned features of P{y) at low T 
one can easily obtain a rough estimate for the superfluid- 
insulator transition point as {J/t)s — 1-05 ± 0.1. How- 
ever, since the precise value is also sensitive to the part 
of P{y) at high fields, j/ > J, a full numerical evaluation 
of the condition (661 is necessary to obtain the exact 
location of the quantum critical point. Using high pre- 
cision data for P{y: T) at low T from Ref. [Ml we find 
{J/t)s ~ 1.00 ±0.01. 

We emphasize an important difference between the 
quantum phase transition we have found here and a stan- 
dard BCS transition. The latter, in the presence of a con- 
stant low energy density of states always yields a finite 
Tc, even though it becomes exponentially small in 1/t for 
small t. In our glassy system the situation is fundamen- 
tally different in that the frustrated interactions suppress 
the density of states around the chemical potential with 
P{y — > 0) — > 0. This quenches the tendency for superflu- 
idity and allows for a superfluid-to- insulator transition at 
a finite value of t, even in the mean field limit of — > oo 
which we consider here. 

This has important consequences for the nature of ex- 
citations and transport properties across the superfluid- 
insulator transition. In particular, the transition to the 
Bose insulator is accompanied by the Anderson localiza- 
tion of lowest energy excitations, whereas higher energy 
excitations remain delocalized relatively far into the insu- 
lator.!^ We believe that the physics revealed by this mean 
field model is relevant for Coulomb frustrated bosonic 
systems which undergo a transition from a superfluid to 
a Bose glass state in finite dimensions. This will be dis- 
cussed in detail elsewhere.!^ 

It is interesting to compare our mean field predictions 
for the phase diagram with the 3D quantum Monte Carlo 
(QMC) simulation results reported in Ref. [TS] The mean 
field predictions for the quantum critical points actually 
match the numerical results surprisingly well. The latter 
were done for the Hamiltonian 

^ = V2)(ri,-l/2) 

(id) 

-t'j2 {bib, +h.c.), 



(67) 



with binary disorder, V-j — ±V' with equal probability. 
Contact with the mean field model ^ is made by re- 
placing the coordination number with N ^ z — 6 for the 
3D cubic lattice, and taking a Gaussian disorder with 
the same variance, V^/z = V'^, as well as a hopping 
h/z = t'. 



Recalling the dictionary ^ , the mean field estimate of 
the superglass to glassy insulator quantum phase transi- 
tion is 



MF 



fV/^z\ (Ul^z 



\ 2t/z 



(68) 



2yi 



4.9 



QMC 



which comes close to the extrapolation of QMC results 
to r = 0. The transition point between superglass and 
non-glassy superfluid is estimated from the static approx- 
imation as 



MF 



2y/z ( y I ~ 2.45 



5,stat 
yl \ QMC 

1j 



3.2. 



(69) 



This indicates that the static approximation overesti- 
mates the stability of the superglass phase, similarly as 
what is known from the mean field version of the trans- 
verse field Ising spin glass. 

The mean field prediction (with static approximation) 
for the interaction-to-hopping ratio {V' /t')T at the tri- 
critical point is rather good, too. 



= 2Vz - 



T,stat 
QMC 



3.55 



3.8. 



(70) 



While the tricritical ordering temperature is overesti- 
mated by a factor of 2 (similarly as in the classical Ising 
spin glass 



MF 



z (T 



2 V t 



T 



\ QMC 



2.2, 



1.1. 



(71) 
(72) 



D. 



Robustness of the phase diagram to random 
field disorder 



In the previous sections we have seen that the model 
([3]) possesses an intermediate phase which is simultane- 
ously superfiuid and glassy. We have determined the 
phase boundaries as instability lines, assuming second 
order phase transitions. Indeed it seems unlikely that 
any of the instabilities could be preempted by a first or- 
der transition. Since the superfiuid to insulator transi- 
tion at {J /t)s is of particular interest, we provide further 
arguments in this section that the parts of the phase dia- 
gram related to the phase boundary of the non-superfluid 
glass remain robust when disorder potentials, i.e. random 
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fields Ci of variance W"^, are restituted to the model. In 
particular we will show that glass and superfluid transi- 
tion lines meet at a tricritical point at finite temperature 
Tt/J and {J/t)T- Further we determine the superfluid 
instability of the glass phase at T = and show that it 
always occurs at a larger ratio ( J/t) than the tricritical 
point, {J/t)s > {J/t)T- This suggests that for any W 
the transition line between non-superfluid and superfluid 
glass is not reentrant as a function of temperature. The 
absence of reentrance in turn suggests that the quantum 
phase transition out of the insulating glass remains sec- 
ond order, independent of the strength of the disorder 
potential. 

The Hamiltonian with a disorder potentials reads 



H 



i<j 



t 

N 



E( 



The disorder potential breaks the Z2 symmetry, therefore 
Qa^b 7^ already in the high temperature phase, where 
it assumes a constant replica symmetric value Qo- The 
glass phase occurs at the Almeida-Thouless instability, 
which is given by: 



/ dy 



Pw{y) 
' cosh'* (/??/) 



1, 



(74) 



where 



Pw{y) 



exp 



2(W2 + J2Q(,) 



and Qa satisfies the self-consistent equation 



Qo = / Fw(2/)tanh^(/3y)dy. 



(75) 



(76) 



The instability towards the superfluid phase is instead 
determined by 



1/1 = j dhPi,,stih) 



tanh(/3/i) 



(77) 



where 



Pi„st(/i) = / dyPwiy) 



cosh /3 ft, 
cosh f5y 



\ 2ho 2 



(78) 



with the Onsager field ho = /3J^(1 - Qo)W^ 

The glass and superfiuid transition lines meet at a tri- 
critical point at {TJJ)t and {J/t)T which are to be eval- 
uated from Eqs. ([^][78|. 

In the limit W/ J ^ 1, one finds the tricritical temper- 



ature (T/J)t 
W/J 



4 J 



and Pxho —J' 3/2, as 



00. 



The superfluid transition at T 
condition 



ms 



dh 



^'inst(^; T 



is given by the 
= 0) 



\h\ 



(79) 



In the hmit W/J » 1, Pi„st(ft; T = 0) is known to have 
a simple structure : 



Pinstih;T = 0) 



v\h\/J^, 







cxp^ 





(80) 



with a smooth crossover between the two limiting forms 
around '-j- = ^^r^^ w + ^iw^)- The value of the con- 
stant 7 = 0(1) can be estimated by the normalization 
condition ^ dhPn^st{h]T = 0) = 1, but will be irrelevant 
below. 

For W/J > 1, {J/t)s and (J/^T both behave as 
vfc '"^wyj'^'' leading order. Their difference scales like 
cJ /W . The coefficient c can be evaluated easily by rescal- 
ing the variables f^Th = h and Pxy ~ y, 



W 

c = lim — \(J/t) 

W/J^oo J 



dh 
h 



pCh:T = 0)-f{h) 



(81) 



where 
and 

f{h)-- 



0) = hm 



VFPinst(ft//?T;T = 0)(82) 



dy smh{h) 



/2tt cosh(y) 



27r 



(83) 



and we have used jSrho — ^ 3/2. We approximate Eq. ( 80 ) 
by extending the formula all the way to /i* and neglecting 
the shift of field h and we get : 



p(ft;r = 0) 



_3_ 



\h\< 



(84) 



Evaluating Eq. (81 1 numerically, using the estimate 
Eq. (84) , one obtains c = 0.231 > 0, establishing that 
{J/t)s > {J/t)T even in the presence of strong disor- 
der. We point out that Eq. (84 1 overestimates Eq. (82 1, 



but this overestimation should be much smaller than 
c = 0.231. 



V. PROPERTIES OF THE SUPER-GLASS 
PHASE 

Having established the phase diagram of the model, 
we now focus on the properties of the bulk of the " super- 
glass" phase. There the interplay between temperature. 
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Fig. 2: Cross-correlation between the local oder parameters 
for superfluidity and glassy order, respectively. The correla- 
tions are evaluated from Eq. ( 85 I in the temperature range 
0.05 < T/t < 0.75 at fixed disorder J/t = 0.724, close to 
the ratio corresponding to the tricritical point. The local 
order parameters are anticorrelated, the maximal anticorrela- 
tion occurs at intermediate temperatures. 



glassy order and superfluid order induce several interest- 
ing phenomena which potentially survive also in finite 
dimensional models of frustrated bosons. In the follow- 
ing, we investigate how the glassy and superfluid orders 
evolve with temperature, and how they are locally corre- 
lated. 



A. Competition between glassy and superfluid 
order 



While the effective transverse field /if = Mt is uniform 
for every site, the frozen longitudinal field, /if depends 
on the site (and on the pure state in which the system 
is frozen). Therefore, the magnetization of the local spin 

Si due to the local field hi = {hf,h^) fiuctuates from 
site to site. It is interesting to study the correlation of 
the local magnetization, whose components are the local 
order parameters of the glassy and the superfluid order, 
respectively. More precisely, we investigate the following 
correlation function: 



C 



M,gEA 



/ dyP{y)m^{y)m'^[y) - Mq_ 



'EA 



Mq- 



'EA 



(85) 



We have evaluated the correlation function ( 85 1 within 



the static 1-step RSB approximation in the center of the 
superglass phase {J/t = 0.724) as a function of tempera- 
ture (0.05 < T/t < 0.75), see Fig.[2] Not surprisingly the 
correlation is negative, since glassy and superfluid orders 
compete with each other. Indeed, one easily checks that 
for every pair of sites it holds that if (sf)^ < (s|) 



then (sf) > (sj). The maximal amplitude of the nor- 
malized correlation Cm,9ea is only of order « 0.1, sug- 
gesting that in the superfluid phase the non-uniformity of 
the two local order parameter fields is actually not very 
strong. It may be that the 1-step approximation under- 
estimates these correlations a bit. The relative weakness 
of the anticorrelations might be the reason why they have 
not been noticed in the quantum Monte Carlo studies of 
Refs.imil 



B. Non-monotonicity of the superfluid order 

In the superglass phase, the glass order parame- 
ter gEA = Qi monotonously decreases with increasing 
temperature, as one should expect. However, surpris- 
ingly, the superfluid order parameter M exhibits non- 
monotonic behavior with a maximum at an intermediate 
crossover temperature T^, as shown in Fig. [3] Below 
r„i, the superfluid order parameter M decreases, anoma- 
lously, when lowering the temperature. Above Tm, M 
decreases with increasing temperature as usual in a stan- 
dard superfluid. 

This phenomenon is related to the anti-correlation be- 
tween glassy and superfluid order discussed in the pre- 
vious section. While on one hand, thermal fluctuations 
tend to diminish both glassy and superfluid order, there 
appears to be a low temperature regime T < T^, where 
quantum fluctuations of the superfluid order are domi- 
nant. Due to the competition between the glassy and the 
superfluid order, the thermally induced decrease of the 
glassy order enhances the superfluid order. This effect 
dominates over the direct thermal effects on the super- 
fluidity. 

It seems natural that it is the superfluid order which 
undergoes such non-monotonic behavior, rather than the 
glassy order. Indeed, we expect the latter to react less 
sensitively to the diminution of quantum fluctuations due 
to decreasing transverse fields 

We note that also the local order parameter correla- 
tions Cm,(jea exhibit a non-monotonous behavior within 
our static 1-step approximation, as shown in Fig. [2j The 
absolute value of CM,qEA increases with temperature at 
very low temperatures, and decreases at higher temper- 
atures. This can be seen again as a consequence of the 
non-monotonicity of the superfiuid order. At fixed T, 
the larger the stronger the normalized anticorrelation 
CM^q^j^. Since = tM initially increases with T, it 
is natural to expect an increasing CM,qEA nntil eventu- 
ally thermal fiuctuations become dominant and diminish 



VI. DISCUSSION 

In this paper we have analyzed a fully connected mean 
field model. The full connectivity is not a real limita- 
tion, however. Indeed, one can generalize the model to a 
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Fig. 3: (color online) The order parameters in the super- 
glass phase as a function of temperature 0.05 < T/t < 0.79 
at the fixed disorder J/t — 0.724. The blue dashed line in- 
dicates the Edwards- Anderson order parameter Qi (in 1-step 
approximation), which monotonously decreases with increas- 
ing temperature. The superfiuid order parameter M (red solid 
line) exhibits non-monotonic behavior. The long time (static) 
on-site charge correlation R, (green dotted line) becomes 1 in 
the disordered high T phase. 



highly connected Cayley tree. While this does not affect 
the thermodynamics of the model, this generalization al- 
lows for the study of localization and delocalization of 
excitations, since this model is endowed with a notion 
of distance. The analysis of localization properties is 
of particular interest in the vicinity of the superglass- 
to-insulating glass quantum phase transition, where the 
boson system collectively delocalizes into a superfiuid at 
low energies. The nature of higher energy excitations in 
the insulator are crucially affected by the suppression of 
low energy states in the glass, leading to a non-trivial ex- 
citation spectrum at the glassy SI transition. The details 
of this analysis will be reported elsewhere.'^ 

What features of the mean field model should be ex- 
pected to carry over to finite dimensions? In the present 
model we find a genuine insulating phase at T = 0, 
which suppresses the superfluidity, due to the strong self- 
generated on-site disorder within the glassy phase. A 
crucial ingredient for the suppression of superfluidity is 
the linear pseudogap within the glass phase. A very sim- 
ilar pseudogap is known to occur in disordered Coulomb 
interacting systems, where it is due to unscreened 1/r in- 
teractions between charged particles. This Coulomb gap 
may well be of importance in strongly disordered super- 
conductors and play a significant role in the competition 
between glassy insulating behavior and superconductiv- 
ity. In particular, in materials with strong negative U 
centers, one may think of preformed electron pairs con- 
stituting hard core bosons which interact with Coulomb 
repulsions.^ The power law suppression of the low energy 
density of states makes it likely that the superfiuid con- 
densate is entirely destroyed once the hopping becomes 



too small. For short range interactions the density of 
states is merely reduced at low energy, but does not tend 
to zero. On a Cayley tree of very large connectivity this 
will always lead to delocalization, unless the hopping t 
is scaled down logarithmically with the connectivity. In 
finite dimensions, however, sufficiently strong disorder is 
known to suppress superfiuidity,'231 and thus one may ex- 
pect that at sufficiently large ratios J/t, the disordered 
boson model will localize due to spontaneously created, 
frozen-in local fields. Such a conclusion may be sugges- 
tive from a straight extrapolation of the quantum Monte 
Carlo results of Ref. [TJ to T = 0, but it seems difficult 
to exclude a scenario in which Tc becomes merely expo- 
nentially small with J/t. A more careful analysis will 
be necessary to settle this question in finite dimensional, 
short range interacting glasses. 

As for the coex istence phase, the " superglass" , the nu- 
merical datcP^'^ provides evidence that it exists also in 
finite dimensions. It would be interesting to confirm 
and quantify the local anticorrelation of order parame- 
ters in such simulations. From our mean field analysis 
one expects that the anticorrelation is in fact relatively 
weak. A further non-trivial prediction with measurable 
consequences is the non-monotonicity of the superfiuid 
order parameter, which should translate into an equiv- 
alent non-monotonicity of the superfiuid stiffness as a 
function of temperature. This non-monotonicity has its 
origin in the softening of the glassy order at low T, a fea- 
ture which may potentially survive in finite dimensions, 
especially when the lattice connectivity is large, or the 
interactions are not too short ranged. We should caution 
though that we obtained this effect by employing a static 
approximation and a replica symmetry breaking at the 
one-step level only. However, we believe that it is a real 
feature of the model. 

As discussed earlier, a number of experiments have al- 
ready shown promising indications of possible coexistence 
of glassy order with superfiuidity. Wc hope that our anal- 
ysis will help to unambiguously identify such phases in 
experiments. Note that finding an experimental system 
exhibiting a glassy superfiuid-insulator transition might 
also be of great interest to study the intricate interplay 
of interactions and disorder with respect to glassy er- 
godicity breaking, and quantum ergodicity breaking, i.e. 
Anderson localization. 
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APPENDIX: Derivation of RS free energy 



With standard replica triclpSl we can get the RS free 
energy: 



— Km — log TrT cxp 



4 



[ f dTdT'R^{T,T') 

Jq Jo 
q(Y1 fdtsl{T)\ 



2 



f dr f dr'R{Ty )sl{T)sl{T') 
a Jo Jq 



+tf3j2 fdrMlslir) 

n Jo 



(86) 



Under static approximation: R{t,t') = R, we have 
/?/ - -■^Q^ + ^M^+'^R^ (87) 



- lim — log TtT cxp 

n->-0 n 



J2/3- 



4 



iR-Q)Y,{f dtsl 

n Jo 



According to Hubbard-Stratonovich transformations, we 
linearize the quadratic terms (^J2a lo '^^'^ai''')^ 
^ dts^^r)^ by introducing extra fields yo £^nd y/j: 



2 2 
2 , M_^2 



tB „ 
— 
2 



(88) 



4^4 

i?yo log / i?2/i^Trexp [/3(yos" + j/fls" + tMs")] 



— T^Q H 7—^ + ^ 



Dyo log / T^yij cosh /3 V (^o + J/fl) + i2M2 



One can get (36 1 following the similar steps above. 
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